GENERAL PURPOSE ANALOG COMPUTATION 



APPLICATION STUDY; 6 



CHEM I C AL / P ETRO -CHEM I C AL 



ANALOG COMPUTER STUDY 
OF A SEMI-BATCH REACTOR 



INTRODUCTION: This paper considers a reactor design application of 
the PACE® TR-48 General Purpose Analog Computer, viz: the chlorina- 
tion of benzene in a semi-batch reactor. The major objective of the 
study will be to relate the maximum production of intermediate chloro- 
benzenes to the chlorine feed rate and the time of reaction. Another 
objective consists of illustrating how experimental kinetic data can be 
fitted on such a computer to obtain reaction rate constants. 

Since the majority of both inorganic and organic chemical products of 
the chemical industry are formed by one or more chemical reactions, 
the chemical engineer is continually faced with the problem of designing 
flow, batch, or semi-batch reactors. In recent years, the cost of design- 
ing and optimizing a reactor has been sharply reduced through the use of 
the general purpose analog computer. The computer is used to simulate 
the behavior of the process in question. Proposed changes in the process 
and its eventual optimization can then be accomplished quickly, and with- 
out expensive equipment modifications. 



THE PHYSICAL SYSTEM 

The chlorination of benzene (CgHg) produces 
monochlorobenzene (CgHgCl) , dichlorobenzene 
(C6H4CI2), and trichlorobenzene (Cg^Cl^) 
through successive, competing chemical reactions: 



C 6 H 6 * C1 2 
C 6 H 5 C1 ♦ Cl 2 

C 6 H 4 C1 2 + C1 2 



-*C 6 H 5 C1 + HC1 (1) 



^►C 6 H 4 C1 2 + HC1 (2) 



-►C 6 H 3 C1 3 + HC1 (3) 



These reactions are carried out in a lead-lined or 
iron vessel, which is shown in Figure 1, with ferric 
chloride (FeC^) as a catalyst. The vessel is fitted 
with cooling coils, since the chlorination reactions 
are exothermic, and a reflux condenser. The pur- 
pose of the reflux condenser is to return vaporized 
chlorobenzenes to the system, while allowing the 
hydrogen chloride and excess chlorine vapor to 
leave the system. In order to maintain the react- 
ing mixture at a uniform temperature and to mini- 
mize mass-transfer effects, the reacting mixture 
is well agitated. 




+> HCI, EXCESS Cl 2 



REFLUX CONDENSER 



DRY Cl 2 FEED 



Figure 1. Semi-Batch Reactor Diagram 

No chlorobenzenes can be formed until some of the 
chlorine gas, which is being bubbled through the 
liquid benzene, has gone into solution. The amount 
of chlorine which is capable of going into solution 
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is limited by the solubility of chlorine in the react- 
ing mixture, which depends on the operating condi- 
tions of the reactor. Since the chlorobenzenes are 
formed successively, the concentrations of mono- 
and dichlorobenzene will each pass through a def- 
inite maximum at some time during the process. 
Eventually, all of the benzene will have reacted to 
pure trichlorobenzene (and other chloro-hydrocar- 
bons, which have been neglected). 

The following restrictions can be made on the system: 

(1) No liquid or vapor hold-up occurs in the 
reflux condenser. 

(2) The system operates under isothermic and 
isobaric conditions. 



N(t) = Moles per unit volume 

k = Rate constant in volume per unit mole 
per unit time. 

and the subscripts represent the following: 

C = Chlorine 

M = Monochlorobenzene 

D = Dichlorobenzene 

T = Trichlorobenzene 

B = Benzene 



(3) Changes in the volume of the reacting mix- 
ture are negligible. 

(4) Hydrogen chloride gas vaporizes and leaves 
the system. 

(5) Negligible mass transfer resistance occurs 
between the gaseous chlorine and the chlo- 
rine in solution (up to solubility limit of 
chlorine). 

THE MATHEMATICAL MODEL 

The equations defining the reaction kinetics of the 
system, given in several references (2,5), are 



dN B (t) 

V —J5 — N B (t) N c (t) 



(4) 



V ^f 1=k . N B (t)N c (t)-k 2 N M (t)N c (t) (5) 
V dJ ^ = k 2 N M (t) N c <t)-k 3 N D (t) NJt) 



(6) 



and 



dN T (t) 
V af— = k 3 N (t) N (t) 



(V) 



1 = Reaction producing monochloroben- 

zene 

2 = Reaction producing dichlorobenzene 
3. = Reaction producing trichlorobenzene 

These equations have the initial conditions: 



N B (0) 


= N o 


(8A) 


N c (0) 


= 0 


(8B) 


N M (0) 


= 0 


(8C) 


N D (0) 


= 0 


(8D) 


N T (0) 


= 0 


(8E) 



Equations 4-7 can be added and integrated to ob- 
tain the benzene material balance equation 

N Q = N B (t) + N M (t) + N D (t) + N T (t) (9) 

If the mass transfer effect of chlorine being trans- 
ferred from the vapor to the liquid phase is ne- 
glected, the chlorine concentration in the system 
is described by the material balance equation 




dt - 



moles of 
chlorine in 
solution 



moles of 
chlorine 
supplied 



x(t) 



moles of 
chlorine 
reacted 



(10) 



where 



V = Volume of reacting mixture where x(t) represents the moles of chlorine reacted 

and r denotes a constant chlorine feed rate in moles 
t = Time per unit time. The chlorine reacted can be obtained 
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(in terms of the chlorobenzene concentrations) from 
the stoichiometry of the system, namely 

x(t) = N M (t) + 2 N D (t) + 3N T (t) (11) 

The solubility of chlorine in the reacting mixture 
was obtained by considering the total pressure above 
the liquid, w 1 which is equal to the sum of the par- 
tial pressures, p, of the inidividual components in 
the vapor phase. Thus 



tt = P c (t)+P B (t)+p M (t) + p D (t) + p T (t) (12) 



Since its tenure in the system is only momentary, 
the partial pressure of hydrogen chloride was ne- 
glected.* By applying Raoult's Law, which states 
that the partial pressure is equal to the product of 

the mole fraction 1 ^ 7 



' N Q + N c (t) 



, and the vapor 



pressure, p , equation 12 becomes 



Nc(t) o N B (t) o | N M (t) o 

" N o +N c (t)Pc + N 0 +N c (t)PB + n o^ n c (+)Pm + 



N D (t) 0 _ N T (t) 0 

N o**" N c (t) P ° N o +N c (t) PT 



(13) 



This equation can be solved for the maximum 
chlorine concentration, NQ(t), to obtain 



N c (t) = 



Np 
P C »-7T 



7T— 



Ne(t) 



Np(t) 
N 0 



(14) 



N T (t) e 



Since the vapor pressures are functions of tem- 
perature only and the system is isothermal and 
isobaric, the maximum chlorine concentration de- 
pends only on the chlorobenzene concentrations. Any 
excess chlorine supplied to the system leaves it 
through the reflux condenser. 

A more convenient form** of the system equations 
may be obtained by defining 



h = 



kN 0 



V 

0 = h t 



(15) 
(16) 



*System operated above critical temperature of HC1. 
**The equations were written in terms of 0, , and because experi- 
mental data was available on the ratios of the rate constants; and k] was 
unknown. 



C(0) = N r l0)/N n 




B(0)= N B (0)/N o 


(18) 


M(0) = N (0)/N o 
M 


(19) 


D(0)= Nd(0 ) /Nq 


(20) 


T(0) = N T (0)/N o 


(21) 


X(0)= x(0)/N O 


(22) 



G =• 



N 0 h, 



(23) 



Using these definitions, equations 4 - 7, 9 and 11 
become 



- ± |f 1 = -B(0)C(fl) 



(24) 



dD(fl) 



dMm =B <e)C(0)-[-j^] M(6)C(0) (25) 

= [y] m{0)c(0)- [-jp] D(e)Cie) (26) 



dT(0) 
d0 



)(0)C(0) 



1= B(0) +M(0) -4-D(0) +T(0) 
X(0) = M(0) + 2D(0) + 3T(0) 

Equation 10 can now be represented as 

C(0) = /G d0 -X(0) 
o 

which can be differentiated to obtain 



dC(0) dX(0) 
= G - 



(27) 
(28) 
(29) 

(30) 



(31) 



d0 d0 

or, by substituting 

1|^.=G-[b(0)C(0)+^) M(0)C(0)+(-^)(32) 
D(0)C(0)] 

Equation 32 is the best form of the chlorine con- 
centration equation for computational purposes. 

The operating temperature and pressure used in 
this investigation are 55 deg. C — which is typical 
of the operating temperature used in industry for 
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this system (4) — and 2 absolute atmospheres — 
which was assumed. Based on these operating 
conditions, the saturation concentration of chlor- 
ine, C(0), can be assumed constant and equal to 
0,120. The proof of this fact is shown in Appendix A. 
This means that the constraint on equation 32 is 



C(6)$ C = 0.120 (33) 



PROGRAMMING THE COMPUTER*** 

This problem, which was solved on a ±10 volt 
PACE® , TR-10 General Purpose Analog Com- 
puter, presents no magnitude scaling problems. 
From the stoichiometry of the system it is obvious 
that B(0), M(0), D(0), and T(0) cannot exceed unity 
and that X(0) cannot exceed three. The chlorine 
concentration, C(0), is limited to 0.120; therefore 
its maximum value is obvious. These maximum 
values and their associated scale factors and com- 
puter voltages are tabulated in Table I. Table II 
lists the systems parametric data. 



TABLE I. MAGNITUDE SCALING 



The scaled voltage equations, which are based on 
the information contained in Tables I and II, are 



J-UtU**\ QOMlkoc] / k 3 \ [|0P][50C] (36) 
dr L I0T J-^J 10 io 



[iop] [50c] 

5k,£ t 10 



[i0b][50b] 
10 



k 2 \ [i0m][50c] /k,\ [lO d] [50 c] 



10 



(- 



10 



(37) 



[?..5X] = 0.250 [l0M] + 0.500 [|0d] + 0.750 [|0d] (38) 



(39) 



Physical 



Estimated Scale Factor Computer 



Variable 


Maximum Value 


Volts Per 


Variable 


(Dimensionless) 


(Dimensionless) 


Dim. Unit 


(Volts) 


B 


1.0 


10 


[10 B] 


M 


1.0 


10 


[io m] 


D 


1.0 


10 


[10 D] 


T 


1.0 


10 


[io t] 


X 


3.0< 4.0 


2.5 


[2.5 X] 


c 


0.12< 0.2 


50 


[50 C] 



TABLE II. PARAMETER DATA 
,(2) 



• = 0.125' 



= 0.00417 



(2) 



Dimensionless 



Dimensionless 



0<G<0.008 (ASSUMED) Dimensionless 



**lt is assumed that the reader is familiar with fundamentals of analog 
computer programming. 



[50 c] < [50 c] =6.00 volts (40) 

where machine time, r, is defined in terms of the 
time scale factor /3, and 6 as 



The computer diagram implementing these equa- 
tions is shown in Figure 2 and the components 
used are summarized in Appendix C. Potentiometer 
and amplifier assignments are tabulated in Tables 
III and IV. The time scale factor chosen for this 
problem was 0.05 (seconds of machine time per 
dimensionless unit). This choice was based on the 
fact that the feedback to the integrators had to be 
attenuated by responsible potentiometer settings in 
order to run the problem in a reasonable length of 
computer time (20 4 t ^ 100 seconds). 

When fitting experimental data to determine rate 
constants, it is advantageous to plot the results in 
time, t, rather than in the dimensionless variable, 
0, since the kinetic data will be in terms of time. 
Time, which is related to 0 by equation 16, can be 
obtained from the equation 



dt_ 
d0 



I 

hi 



k, N 0 



(42) 
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+10 



Ol 




Figure 2. Computer Diagram 



TABLE HI 



ELECTRONIC ASSOCIATES INC. 
EDUCATION a TRAINING GROUP 

TR-IO POTENTIOMETER ASSIGNMENT SHEET 
PROBLEM ^^^'-RsfeSL jcW g&g^. 
DATE 



POT 

NO. 


PARAMETER 
DESCRIPTION 


SETTING 
STATIC 
CHECK 


STATIC 
CHECK 
OUTPUT 
VOLTAGE 


oil 1 1 11Mb 

RUN 
NUMBER 1 


NOTES 


POT 
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1 


2 


i/fib/9 




o. 4-oo 







2 


3 










5* 


3 


4 










> 


4 








5 




o. £"Oq 








5 








6 






v » t > s g 






6 


7 


v. o fvj ^> \ r*— ' > 








^ 


7 


8 




on sro 








8 






- > 


9 




o . ZS"D 









9 








10 




CM OO 


O. OOO 




> 


10 


II 




(D , too 








I 1 






^ 


12 


L< M (T 








> ftPPf^iMATe 


12 








13 












13 


14 


L i m rr 








> APPROXIMATE 


14 








15 












15 


16 












16 


17 












17 


16 












18 


19 












19 


20 












20 


21 












21 
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23 


FS1 


LAP 
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23 


24 




UP 


up 
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24 



M673 
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TABLE IV 



ELECTRONIC ASSOCIATES INC 
EDUCATION a TRAINING GROUP 

TR-IO AMPLIFIER ASSIGN MEN SHEET 
DATE 



AMP 
NO. 


FB 


OUTPUT 
VARIABLE 


STATIC CHECK 


NOTES 


CALCULATED 


MEASURED 


INTEGRATOR 
INPUT SUU 


OUTPUT 


INTEGRATOR 
WPUT SUU 


OUTPUT 


1 


I 






\G, (OO 








2 


I 


- loB 




-iO.OO 








3 








— (o.OO 








4 


s 






10.00 
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i oM 
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TABLE V. TIME OF MAXIMUM PRODUCTION 
VERSUS CHLORINE FLOWRATE 
DATA 



6 Max 
Dimensionless 



G 

Dimensionless 



l/G 

Dimensionless 



M 


D 






22 


264 




0 


160 


430 


0.008 


125 


210 


520 


0.006 


167 


305 


680 


0.004 


250 


585 


1230 


0.002 


500 



whose scaled voltage representation is 
d [lOt/t L ] 



(JT 




(43) 



The symbol t L is the maximum time required for 
the data analysis, and depends on the data being 
fitted. 

A dynamic check of the computer circuit was made 
after completion of successful static check. The 
dynamic check consisted of checking the maximum 
concentration of monochlorobenzene and the time 
at which the maximum occurred, against the com- 
puter results for a special case of reactor opera- 
tion. The special case (FSI -DN, FS2-UP) assumes 
that the chlorine flow rate is infinite and that the 
reaction mixture is always saturated with chlorine. 
Calculations made on this basis indicate the maxi- 
mum value of M(0) will be 0.74 when 0 is 22 (Refer 
to Appendix B). 

RESULTS 

In addition to the dynamic check run, as shown in 
Figure 3, four additional runs were made at values 
of G ranging from 0.002 to 0.008. These results, 
which are plots of B(0), M(0), D(0), T(0), X(0) and 
C(0) versus 0, are shown in Figures 4-7. From 
this information, it was possible to obtain the 
values of 0, 0™ QV at which the maximum amount of 



the intermediate products was present in the sys- 
tem. This data and the corresponding values of the 
flow rate variable, G, are tabulated in Table V. A 
plot of log 0 max versus l/G is shown in Figure 8 
for each of the intermediate products. 

Figure 9, which is a plot of B(t), M(t), D(t), and 
T(t) versus time for three values of k3/k^, serves 
to illustrate how experimental data can be fitted 
using the analog computer. By varying the reaction 
rate constants, which are unknown, the computer 
solution can be made to coincide with the experi- 
mental results to obtain these constants. 
When this trial and error process is performed in 
real time, using an X-Y plotter, it is quite time- 
consuming. However, when high speed repetitive 
operation, which yields up to 60 solutions per sec- 
ond, is used in conjunction with an oscilloscope, a 
transparent overlay of the experimental data can be 
fitted in a matter of minutes, due to the fact that 
the effect of parameter variations on computer 
results can be observed instantly on the high per- 
sistence scope. 

CONCLUSIONS 

From Figure 3-7 it is obvious that the maximum 
reduced concentration of the intermediates is inde- 
pendent of the chlorine flow rate. These maximums, 
which occur at 0.74 and 0.88, can only be altered 
by changes in the reaction rate constants. This can 
be verified by consulting Figure 9. The chlorine 
feed rate, which controls the amount of chlorine 
available for reaction, affects the time at which 
the maximum concentrations occur. As shown in 
Figure 8, when the chlorine flow rate is infinite 
(l/G =0) and the reaction mixture is always satur- 



ated with chlorine, 0 



max 



has its minimum values 



of 22 and 264. When the chlorine feed rate is re- 
duced, # max increases. Obviously, if the feed rate 

of chlorine is zero (l/G-^oo) 0 must be infinite. 

' max 

The curve fitting technique mentioned in the results 
section, which is for the most part qualitative, 
serves only as an illustration of what can and has 
been done on the general purpose analog computer 
in regard to the analysis of kinetics data. 
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Dimensionless Time, 0 

Figure 3. Reduced Concentrations Versus Dimensionless Time For Infinite Chlorine Feed Rate 




O 200 400 000 000 1000 1 200 1400 13O0 I0OO 2000 2200 



Dimensionless Time, 0 

Figure 4. Reduced Concentrations Versus Dimensionless Time (G = 2.0x10"^) 
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Figure 7. Reduced Concentration Versus Dimensionless Time (G = 8. 0x10" ) 




Reciprocal Chlorine Flow Rates, G" 1 



Figure 8. Dimensionless Time of Maximum Mono-and Dichlorobenzene 
Production Versus Reciprocal Chlorine Flow Rate 
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0 10 20 30 40 50 60 



Time, Minutes 



Figure 9. Reduced Concentrations Versus Time, Comparison of Experimental and Computer Results 



APPENDIX A 
SOLUBILITY CONSTRAINT CALCULATIONS 



The vapor pressure data at the 55 deg. C operating temperature, necessary to solve the equation 
I 



C( t) = 



0-B(t)p B ° -M(t)pJ -D(t) p° D -T(t)pf ] 



is tabulated in Table VL 



(A-l) 



TABLE VL VAPOR PRESSURE DATA AT 55 C 



Component 
Benzene 

Monochlorobenzene 
Dichorobenzene 
Trichlorobenzene 
Chlorine 



Vapor Pressure 

327 mm of HG 

48.8 mm of HG 
11,5 mm of HG 
<5 mm of HG 

16.9 ATM 



Data 
Source 

3 
3 
3 
3 
1 



If this data is substituted into equation A-l, it can be seen that for all practical purposes since the total 
pressure is two atmospheres (see Table II), and the concentration variables cannot exceed unity. Based 
on equation A- 2 the range of C(t) is 

I 



C(t) 



PC" 



O.I05< C( t)< 0.134 



(A-2) 



(A-3) 



which is very narrow; therefore, variations in C(t) will have little effect on the results of the problem. The 
average value of C(t), 0.120, was considered sufficiently accurate for this study. 
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APPENDIX B 
DYNAMIC CHECK CALCULATIONS 



If the reacting mixture is always saturated with chlorine, the equations describing the behavior of 
chlorobenzene and benzene are 

dB 

^(f)=-B(f) 



and 



where 



These equations have the initial conditions 



« «,..(«- 4«- ma 



e = ce 



B(0)= I 
and M(0)= 0 



The solution to equation B-l is 



B {£) = e" f 



which can be substituted into equation B-2 to obtain 

d_Mii + £l M , n= .-« 

whose solution is 

M(£) = 



The maximum concentration of M(£) occurs when 



.if 

dM(Q =0= / _ J^2. e kl 



where i; ' is the value of the argument at which the maximum occurs . 



e = -, In 

I -ILL k 2 



k 



This may be substituted back into equation B-8 to obtain the maximum value of M(£ ) 

kg/k, 



k 2 

Since — is one-eighth (see Table II) and C is 0.12 (see Appendix A) 
1 M(€) =0.74 

and d ' = -4- = 22 

C 
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APPENDIX C 

EQUIPMENT REQUIRED TO PERFORM SIMULATION 



The equipment required to perform this simulation was: 

1 - X-Y Plotter 

3 - Quarter Square Multipliers 
17 - Operational Amplifiers 

6 - Integrator Networks 
13 - Potentiometers 

and 

2 - Function Switches 
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NOMENCLATURE 

B = Reduced Benzene concentration, dimension- 
less 

C = Reduced Chlorine concentration, dimension- 
less 

C = Reduced Maximum Chlorine concentration, 
dimensionless 

D = Reduced Dichlorobenzene concentration, di- 
mensionless 

G = Modified chlorine flowrate, dimensionless 

M = Reduced Monochlorobenzene concentration, 
dimensionless 

N = Concentration variable, moles 

Nq = Maximum chlorine concentration, moles 

T = Reduced trichlorobenzene concentration, di- 
mensionless 

V = Volume of reacting mixture, cubic feet 

X = Reduced reacted chlorine variable, dimen- 
sionless 

h = Modified rate constant, seconds 1 

k = Reaction rate constant, cubic feet/mole sec- 
ond 



p = Partial pressure, atmospheres 

pO = Vapor pressure, atmospheres 

r = Chlorine feed rate, moles/ second 

t = Time, seconds 

t L = Maximum time, seconds 

x = Moles of chlorine reacted 

t = Absolute pressure, atmospheres 

0 = Dimensionless time variable, dimensionless 

T = Machine time, seconds 

P = Time scale factor, seconds/dimensionless 
unit 

£ = Dimensionless time variable, dimensionless 
SUBSCRIPTS 

B = Benzene 

C = Chlorine 

D = Dichlorobenzene 

M = Monochlorobenzene 

0 = Initial Value 

T = Trichlorobenzene 

1 = Reaction producing Monochlorobenzene 

2 = Reaction producing Dichlorobenzene 

3 = Reaction producing Trichlorobenzene 
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